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thereby prohibit the existence of propagating vector fields. The induced gravito- electromagnetic Lorentz 
force on a test particle is evaluated in terms of these fields together with the torque on a small gyroscope. 
It is concluded that the analogy of perturbative gravity to Maxwell's description of electromagnetism can 
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1. Introduction 

Einstein's theory of gravitation remains a pinnacle in the evolution of theoretical physics. It offers an 
overarching description of phenomena ranging from the familiar behaviour of Newtonian gravitation to exotic 
astrophysical events at the extremes of space and time. Although its modern formulation is in terms of tensor 
fields on a manifold with a spacetime structure, its physical interpretation often benefits from a choice of 
observer and an appropriate reference frame. One of the traditional methods for extracting information from 
Einstein's gravitational field equations is to exploit the properties of observers in some fiducial background 
spacetime in which gravitational physics is either absent or familiar. This approach has led to various 
approximation in such backgrounds. Further reduction is afforded by a "3+1" decomposition in which 
spacetime tensors are expressed in terms of a field of frames adapted to some local foliation of spacetime by 
spacelike hypersurfaces. More generally frames are afforded by timelike vector fields, the integral curves of 
which model ideal observers. In such a manner it becomes possible to contemplate different limits in which 
matter moves slowly or the gravitational field is weak relative to such observers. It is known that Newtonian 
gravitation follows from such a limit. Within the framework of weak gravity non-Newtonian gravitational 
effects may arise and a number of experiments have been devised in order to detect such phenomena as the 
"Lense-Thirring effect" due to "frame-dragging" produced by the earth's rotation p2| , Q|, ||. The 
nature of this effect may be detectable by a small orbiting gyroscope and is analogous to that produced by 
the torque on a small magnetic dipole in the presence of the magnetic field of a fixed magnetic dipole. Indeed 
the component of weak gravity (additional to the dominant Newtonian gravitational field) responsible for 
this effect is now referred to as the gravito-magnetic field. The sensitivity of recently developed rotation 
sensors may also be increased to detect post- Newtonian effects in the future Jl9{ . 

Several authors have noticed that a subset of the Einstein 

equations when perturbed about flat spacetime can be written in a form that looks remarkably similar to 
Maxwell's equations with the Newtonian gravitational field corresponding to the gravito- electric field and 
mass-currents playing the role of electric currents. Since the laws of electromagnetism are well studied 
and understood this analogy has proved quite fruitful in the gravito- electromagnetic context particularly in 
astrophysical applications. Extended "astrophysical jet-structures" are now thought to have their origin in 
gravito- electromagnetic forces. In [ |l8| the details of astrophysical lensing have been explored in terms of 
parameters in the NUT metric. 

It is also amusing to recall that one of the first theories of post-Newtonian gravitation was formulated 
by Heaviside in direct analogy with the then recently formulated theory of electromagnetism by Maxwell. 
In the language of the Poincare isometry group it predicted that gravitation like electromagnetism was 
mediated by an independent vector field rather than a second degree tensor field associated with the metric of 
spacetime. This difference of course must imply that any analogy between weak gravity and electromagnetism 
is incomplete and most derivations of the gravito- electromagnetic field equations take care to point this 
out. However in our view the caveats are themselves often incomplete and a close examination of various 
derivations of the gravito- electromagnetic equations display significant differences in detail. The difficulty 
in making objective comparisons often arises due to the implicit use of a particular coordinate system 
(usually adapted to a flat spacetime background) or a partial gauge fixing. Indeed the question of the gauge 
transformations induced on the gravito- electromagnetic fields from the underlying gauge covariance of the 
perturbative Einstein equations is usually dealt with rather cursorily. This leads one to contemplate the most 
useful way to define the gravito- electromagnetic fields in terms of the perturbed components of the spacetime 
metric. Different choices are often responsible for the location of odd factors of 4 that permeate the gravito- 
electromagnetic equations compared with the Maxwell equations. Such choices also have implications for the 
form of the induced gravito- electromagnetic Lorentz force (and torque) in terms of the gravito- electromagnetic 
fields that enter into the equation for the motion of a massive point (spinning) particle. In fT(| gravito-electro- 
magnetic gauge transformations are discussed from a perspective different from the one presented in this 
paper. Here such transformations are explicitly related to the gauge symmetry of perturbative gravitation 
and the definitions of gravito-electromagnetic fields in turn induce the notions of gravito-magnetic and 
gravito-electric coupling strengths. 

In this article a tensor description of perturbative Einsteinian gravity about an arbitrary background 
spacetime is first constructed. By analogy with the covariant laws of electromagnetism in spacetime gravito- 
electromagnetic potentials and fields arc then defined to emulate electromagnetic gauge transformations under 
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substitutions belonging to the gauge symmetry group of perturbative gravitation. These definitions have the 
advantage that on a flat background, with the aid of a covariantly constant timelike vector field, a subset of the 
linearised gravitational field equations can be written in a form that is fully analogous to Maxwell's equations 
(without awkward factors of 4 and extraneous tensor fields. It is shown how the remaining equations in the 
perturbed gravitational system restrict the time dependence of solutions to these equations and thereby 
prohibit the existence of propagating vector fields. The induced gravito- electromagnetic Lorentz force on 
a test particle is evaluated by geodesic perturbation in terms of these fields together with the torque on 
a small gyroscope. It is concluded that the analogy of perturbative gravity to Maxwell's description of 
electromagnetism can be valuable for (quasi-) stationary gravitational phenomena but that the analogy has 
its limitations. It has been argued that such limitations are absent in the approach to gravito- electromagnetic 
based on properties of the conformal tensor in a spacetime determined by Einstein's equations. Although 
this reformulation makes no reference to perturbative methods the analogy with the structure of Maxwell's 
equations is less direct. A tensorial description of this formulation is given in Appendix |d|. Throughout 
this article the language of tensor fields as multi-linear maps on vector and co-vector fields is adopted. Co- 
vector fields are manipulated using the exterior calculus of differential forms and Hodge maps. Manifolds are 
assumed smooth and tensor fields sufficiently differentiable as required. Notations based on the tools used are 
summarised in Appendix ^ and some technical computational details are relegated to Appendices |b] and 0. 
In order to facilitate comparisons with other authors certain field redefinitions are discussed in section |7j 
together with the changes induced by them in the gravito- electromagnetic field equations. These alternatives 
are discussed in the concluding section where the salient features of this paper are summarised. 

2. The Maxwell Equations 

Maxwell's equations for the electromagnetic field can be written concisely (in a general spacetime with 
metric tensor g) in terms of the Faraday 2-form F and current 1-form J on spacetime: 

dF = 0, (2.1a) 
SF = J. (2.1b) 



Equation ( 2.1a ) implies that in a regular domain 

F = dA (2.2) 

for some 1-form potential A, and F clearly remains invariant under electromagnetic gauge transformations 
of A which take the form 

A^A + dX, (2.3) 



where A is an arbitrary smooth function on spacetime. Equation (2.1b) implies that the current is conserved 

SJ = 0. (2.4) 

In terms of A, ( |2.1a ) may be written 



AA + d5A = -J, (2.5) 
where A = — (<5d + dS). In the Lorenz gauge defined by the condition 

SA = 0, (2.6) 

this reduces to the Helmholtz wave equation 

AA = -J. (2.7) 

The equation of motion for a (spinless) test particle with mass m and charge q moving along a curve 
C(r), parameterized by proper-time r, with tangent vector C"(r), in an arbitrary gravitational field and 
electromagnetic field F is 

mV C 'C' + ?(ic'F)"=0, (2.8) 
where V is the Levi-Civita connection associated with g. The Lorentz force Fl is defined by 

¥ L = -q(i c ,F) S . (2.9) 
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One can decompose F with respect to some unit-normalized timelikejj] vector field V as 

F = F b Ae + #b, (2.10) 

where iy-e = and iy-^b = 0. Similarly one may decompose J as 

J = pV b +i (2.11) 

where iy-j = 0. Ideal observers may be associated with the integral curves of V . 

The Maxwell field equations ( 2.1a ) and ( 2.11 ) can then be written in (3+l)-form in terms of e, b and V 

as 

dy#b + n v A e = 0, (2.12a) 
Dve + £v#b = 0, (2.12b) 
dy#e-nvAb = /j#l, (2.12c) 

Dy-b - £ y #e = #j, (2.12d) 

where the notation is given in the Appendix |a|. Such equations exhibit possible "pseudo-sources" measured 
by non-inertial observers associated with non-parallel vector fields V. 

When V is a parallel vector field (e.g. associated with an observer defining an inertial frame in flat 
spacetime) these reduce to the familiar Maxwell equations relating electromagnetic fields to their sources 
written in terms of differential forms Q. The metric tensor permits one to define the vector fields: E = e", 
B = b", J = j", and (for V parallel) the Maxwell equations take their more familiar form 

divB = 0, (2.13a) 

curlE+— =0, (2.13b) 

divE = p, (2.13c) 
<9E 

curlB- — =J. (2.13d) 
Similarly decomposing the 4-velocity C with respect to some (unit timelike parallel) vector field V as 

C'= . = =(^ + v), (2.14) 

v 1 - g( v > v ) 

and assuming that the magnitude of v is small {v 2 = g(v, v) -C 1), 

C = V + v + 0(v 2 ), (2.15) 
the equation of motion for a charged non-relativistic particle becomes 

^ = -(E + vxB). (2.16) 
dr m 

3. Perturbative Gravitation 

The analogy between gravitation and electromagnetism to be discussed follows from a perturbation of 
the gravitational field about some "background" spacetime geometry. We explore the constraints on this 
geometry in order to execute this analogy to its fullest extent. 



-^uch a vector satisfies g(V, V) = 



— 1 with the choice of signature for g chosen throughout this article. 
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3.1. The Perturbed Einstein Equations. 

A perturbative approach to Einstein's theory of gravitation can be based on a formal linearization of 
spacetime geometry about that determined by some solution of Einstein's equations for the spacetime metric 
tensor. A generic perturbation will be identified with a class of linearizations that belong to the tangent space 
to the space of Einstein solutions. Following |2(J it is convenient to introduce a 5-dimensional manifold M 
that can be foliated by hypersurfaces belonging to a one parameter (e) family of spacetimes. The geometry 
of each leaf M e of the foliation is determined by some second degree symmetric tensor field g(e) on M. where 
e E [—1,1] that restricts to a Lorentzian metric tensor field on each leaf. Furthermore it is asserted that all 
points of one leaf can be connected to all points on a neighbouring leaf by a one parameter diffcomorphism 
ip e : M. — > M.. Thus a tensor field T(0) on Mq can be related to a tensor field T(e) on M c according to 

^_ e T(e) = T(0) + eT v + O(e 2 ) (3.1) 

where, for some nowhere vanishing vector field V that is nowhere tangent to the spacetime leaves in M, 

<p € = exp(eV) (3.2) 

induces ip e on tensors and 

TV = C v T(e)\ e=0 . (3.3) 
In a local chart with coordinates e} adapted to the foliation one may write: 

The tensor Xy is said to be a linearization of T(e) about T(0) with respect to a choice of the vector field V. 
Since the leaves M e are diffeomorphic it is natural to identify points on distinct leaves that lie on the same 
integral curve of V. Different choices of V correspond to different identifications. If Ty t and Ty 2 are distinct 
linearizations then by construction their difference is generated by a vector field X on Mq: 

T Vl ^ T V2 = f Vl + C x T(0). (3.5) 

This may be called a gauge transformation of T\> 1 induced by X. If CxT(0) = for all X then T(e) is gauge 
invariant. In general there is no preferred positive (or negative) definite metric on Mq that enables one to 
assign a natural norm to quantities that are not gauge invariant in this sense. 

Linearized gravity proceeds by writing the covariant physical metric tensor g(e) field so that: 

cp^g(e) = ,g(0) + eg v + 0(e 2 ) (3.6) 

By common abuse of notation this is simply written; 

g = g + h (3.7) 

where h is of order e and higher order terms in e are subsequently neglected. This abbreviated notation 
hides the choice of V in the definition of h and an alternative choice arises from the gauge transformation 

h^h + eC x g (3.8) 

for any vector field X on the spacetime with background metric g. Similarly we introduce the abbreviated 
notation f = ifi- e T(e), T = T(0), f = ef v , for some tensor T (thus t is of order e). 

The contravariant physical metric tensor field G can be similarly written in terms of the induced con- 
travariant background metric tensor G and a contravariant perturbation tensor H 

G = G + H. (3.9) 

In a smooth local basis of vector fields on spacetime {A" a }, with dual cobasis {e b } such that e b (X a ) = S b a 
for a — 0, 1,2,3, the induced perturbation tensor H can be written as 

H=-h ab X a ®X b (3.10) 

with g ac g cb — S b a (where g ab are the components of G in the above cobasis) and X a = g ab X},. The metric- 
dual of vector and 1-form fields, and the associated raising and lowering of indices are defined with respect 
to the background metric tensor. Likewise any operations that depend on the metric tensor (e.g. the Hodge 
map) will be defined with respect to the background metric unless indicated otherwise. 
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To derive the perturbed Einstein tensor in terms of h (in some gauge) write the Levi-Civita connection 
V with respect to the physical metric g in terms of the Levi-Civita connection V with respect to g so that 
for any vector field X on spacetime: 

V x =V x +~fx, (3-11) 

Since V and V are torsion free 

lx Y = lY X, (3.12) 
for arbitrary vector fields X and Y. It follows that 

7a/ = 0, (3-13) 
for any function / on spacetime. It is convenient to define a tensor 7 by 

1 {X,Y,a) = a{ lx Y), (3.14) 

for any 1-form a, then 

lx Y = 1 {X,Y-) (3.15) 

and 

lx a=- 1 (X,-,a). (3.16) 
Since V and V are compatible with respect to g and g respectively 7 can be written in terms of h as 

l(X,Y,a) = ^{(Vx/i)(l>S) + (V Y h)(X,J) - (V at h)(X,Y)}, (3.17) 
for arbitrary X, Y and a. 

The curvature operators R and R for the connections V and V respectively are defined in the usual 
manner as 

Rxy = [Vx.Vy]- V [X ,Y], (3.18a) 
Rxy = [Vx,Vy]-V [x ,y]. (3.18b) 

Thus 

Rxr = Rxr + Rxy, (3.19) 
and it follows that the perturbed curvature operator R is given by 

Rx y = Vx7y - Vy7 x +7 x Vy- 7yV x - 7[a,y] ■ (3-20) 
When acting on any vector field Z this simplifies to 

RxyZ=(Vx1)(Y,Z,-)-(V y7 )(X,Z,-), (3.21) 
and the perturbed curvature tensor R is defined by 

R(X, Y, Z, a) = a(K x Y Z) = (V x7 )(r, Z, a) - (Vy 7 )(X, Z, a). (3.22) 

The perturbed Ricci tensor Ric follows by contraction 

Ric(X, Y) = R(X a , X, Y, e a ). (3.23) 
It is useful to introduce the trace-reverse map \x on covariant degree 2 tensors so that 

H(T) = T - \ Tr(T) g (3.24) 
where Tr(T) = T{X a , X a ) in any basis. The trace-reverse of h is denoted by ip: 

i> = fj,(h). (3.25) 

After some calculation Ric can be written in terms of ip and the Laplacian operator Lap = V-V as 

Ric = -i^(Lap V) + SymVV--0 - + S^, (3.26) 
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where 

Ctp(X,Y) = R(X a ,X,Y^ a ), (3.27a) 



1 
2 

and the convenient notation 



S^X, Y) = -(Ric(4, Y) + Ric(V^, X)) (3.27b) 



i>x=i>{X,-), (3.28a) 

r =ijj X a =ip{X a ,-) (3.28b) 

is used. 

Note that C^, and have the same trace, namely 

Tr(C^) = Tr(S^) = Ric(V|, X a ). (3.29) 
The perturbed curvature scalar 1Z follows as 

Tl = Tr(Ric - S t + ^Kip), (3.30) 

and with Ein = Ric — \gfl = Ein + Ein the perturbed Einstein tensor Ein is 

Ein = Ric - -Kh - ^Tlg, (3.31) 

or 

Ein = -iLap^ + Ai(SymVV» - fi(C^) + S v , - ^Tltp (3.32) 

in terms of ip. 

One may now express the perturbed Einstein equation in terms of ip by writing the physical stress energy- 
momentum tensor T as 

t = r + t, (3.33) 

where the background stress energy-momentum tensor T acts as a source for the background metric via the 
background Einstein equation 

Ein = kT, (3.34) 
and k — 8ttG in units where c = 1. The perturbed Einstein equation 

Ein = kT, (3.35) 

becomes 

+ fi(SymVV^) - n(C t )+S^ - ^Kip = nt. (3.36) 
In the next section this equation is simplified by exploiting a gauge symmetry of the Einstein equations. 

3.2. Gauge Transformations and the Transverse Gauge Condition. 

A gauge transform of h has been defined as a substitution of the form 

h^h + C v g, (3.37) 

where Cy is the Lie derivative with respect to some vector field V that maintains h perturbative with respect 
to g. This substitution is used to determine the induced gauge transformation of tensors defined in terms of 
h (keeping the background geometry fixed) . Thus the induced gauge transformation of f x follows as 

1 X ^1 X +V X (V), (3.38) 

where, for any vector field X the tensor derivation T>x{V) is defined by 

V X (V) = [C v ,Vx] -V Cv x. (3.39) 

The operator T>x{V) provides a useful tool when performing calculations involving both Lie derivatives and 
covariant derivatives (Ej). 
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It follows that the perturbed curvature operator transforms as 

K X Y^RxY+n XY {V), (3.40) 
where the operator 1Zxy(V) is defined by 

Ttx y(V) = [Cv , R-x y] ~ R-c v x y — Rjc c v y- (3-41) 
A contraction shows that the perturbed curvature tensor transforms as 

R^R + CvR, (3.42) 

from which one deduces 

Ric i-> Ric + £ y Ric (3.43) 

and 

K ^ K + C V U. (3.44) 
It follows that the perturbed Einstein tensor exhibits the induced gauge transformation 

Ein i-> Ein + £ y Ein. (3.45) 

This behaviour of Ein dictates the behaviour of any (phenomenological) stress energy-momentum tensor 
under an induced gauge transformation in order to maintain the gauge covariance of (3.35). 
It can be similarly shown that if ip i— > ip = ip + ix{Cyg) the divergence of ip transforms as 

V-ip>-> V> = V> + (Lap Vf + Ric(K -). (3.46) 

Thus, if for some ip, one chooses V to satisfy the differential equation 

(LapV) b +Ric(F,-) = -V-ip, (3.47) 

then this imposes the gauge condition 

V-^ = 0, (3.48) 

referred to as the transverse gauge. Further gauge transformations, generated by additional vector fields W 
satisfying the linear differential equation 

(Lap Wf + Kic(W, -) = 0, (3.49) 

maintain ip divergenceless. 

In the transverse gauge the perturbed Einstein equation immediately simplifies to 

~Lapil> - H{C^,) - -Hip = k±. (3.50) 

It should be noted that this equation holds in any (background) metric that satisfies the (background) 
Einstein equation ( 3.34 ). 

In the case where the background metric is Ricci-flat (or has a cosmological constant), the perturbed 
Einstein equation simplifies to 

-^LapV-Cv, = nf (3.51) 

which may be compared with the component form p"7j ] and the abstract index form p5| . 
If the background metric is flat this simplifies further to 

Lap-0 = -2Kr. (3.52) 

This equation is analogous to the Helmholtz equation for the electromagnetic potential 1-form in the elec- 
tromagnetic Lorenz gauge and leads to the prediction of propagating gravitational perturbations. 
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4. Gravito-Electromagnetism 

The perturbed Einstein equation about a flat background spacetime metric has a form that is similar to 
the Helmholtz equation for the electromagnetic potential. However the perturbed gravitational potential ip 
is a symmetric type (2,0) tensor field, whereas the electromagnetic potential A is a type (1,0) tensor field 
(a 1-form). The description of Maxwell's equations in terms of electric and magnetic fields suggests that 
ip be decomposed relative to some fiducial vector field and split into a frame dependent 1-form (analogous 
to the 1-form potential in electromagnetism) plus an extra non-Maxwellian tensor field. We motivate this 
definition by inducing a gauge transformation described above on this 1-form and comparing the result with 
the structure of the electromagnetic gauge transformation (Section ||) of the 1-form potential. This will also 
show that restrictions must be placed on the fiducial vector field if a close analogy with electromagnetism is 
to hold. 

4.1. Gravito-Electromagnetic Gauge Transformations. 

The trace-reversed perturbation ip can be written, relative to some unit-normalized (to zero order in e) 
timelike vector field £, (g(£, £) = — 1 + 0(e)) and general background metric tensor g as 

V = ^5-^®^-e b ®^-S € . (4.1) 

The tensor ip has been split into a 1-form part 

1> e = (4.2a) 



with £ component 

and a spacelike tensor part satisfying 
Since 



= ^(0=^.0. ( 4 - 2b ) 
£ e (£,-)=0. (4.2c) 



iP^iP + ^Cyg) (4.3) 
under a gauge transformation h i— > h + Cyg generated by V, a contraction with £ gives 

rl>t~1>t + V)) + iyd£ b - V-V? + (C^vf. (4.4) 
With a special choice for V of the form 

V = -A£, (4.5) 

where A is any smooth function of order e on spacetime, the induced transformation on the 1-form ip^ becomes 

^^^+dA + A(6 e £ b -^), (4.6) 

where 0£ is the expansion of £ and A.^ is its acceleration (Appendix [A|). If A(9^£ b — «4je) is of higher order 
in e than A the last term can be neglected and the gauge transformation of ip£ simplifies to 

ipt: i-> tp£ + dA. (4.7) 

This suggests that ip£ be interpreted as the analogue of the 1-form potential A in electromagnetism and 
henceforth ip£ will be ref erre d to as the gravito- electromagnetic l-form potential. 
Similarly, projecting (4.5) with IIj (Appendix |a|) one finds 

S £ ~ S f - U^C v g + 2(V- V + g(£, V 5 F))g 5 . (4.8) 



Hence with V as in (4.5), transforms as 

E € ~££+2AK-fe^), (4.9) 

where er^ is the shear of £ (Appendix ^). If A(cr £ — |0£g£) is of higher order in e than S £ the last term is 
negligible and is gauge invariant under such transformations. 



The class of gauge transformations generated by — A£ such that ip£ transforms as in (4.7) with invariant 
will be said to contain gravito- electromagnetic gauge transformations. Clearly if £ is parallel (i.e. V£ = 0), 
or parallel to zero order in e (since A is of order e), such transformations are guaranteed to be gravito-electro- 
magnetic gauge transformations. For a given h one may therefore define a class of vector fields {£}, members 
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of which are equivalent if they can be used to generate gravito-electromagnetic gauge transformations. Such 
fields will be said to define gravito-electromagnetic frames of reference. Not all background spacetimes will 
permit the existence of such vector fields (e.g. the black- hole Schwarzschild spacetime), however they are 
guaranteed to exist in Minkowski spacetime (and interestingly also in the Einstein static universe). In the 
next section these vector fields will be used to define a class of gauge equivalent gravito-electromagnetic fields. 

4.2. The Gravito-Electromagnetic Field Equations. 

On a flat background the linearized Einstein equation can be written as 

LapV = -2KT, (4.10) 

where the transverse gauge condition 

V-V> = 0, (4.11) 



has been imposed. Contracting ( 4.10 ) with a unit-normalized timelike vector field £ that can be used to 
define gravito-electromagnetic gauge transformations gives 

A^ = 2nJt, (4.12) 

where the mass-current J7j is defined by 

Jt = -tfo-). (4.13) 



Contracting (4.12) on £ gives 

A0 C = -2 m , (4.14) 

where the mass-density is defined by 

Pi = -jm=T{u)- (4.i5) 



Similarly acting on ( 4.10 ) with the projection operator 11^ and using (4.1) 

A0 c g 5 - Lap S c = -2kT 6 . (4.16) 



With the aid of ( 4.14 ) this becomes 

LapS 4 = 2/e(7' c -p c g c ), (4.17) 
where is the spacelike (with respect to £), part of the stress energy- momentum tensor 

±z = n ? r. (4.18) 

The gravito-electromagnetic analogue of the Faraday 2-form is now defined as 

Ft = d^, (4.19) 

hence 

= 0. (4.20) 

The tensor J-^ is invariant under transformations in the class of gravito-electromagnetic gauge transforma- 
tion defined above. The gravito- electric field and the gravito-magnetic field 3^ follow from a (3+l)-split 
(with respect to £): 

^=^A£ ? + #B e . (4.21) 
With written in terms of its timelike and spacelike parts: 

^ = -^ b + *e, (4-22) 

where if*£ = 0, one can write the gravito-electric field as 

= d£<f>£ - (4.23a) 

and the gravito-magnetic field as 

B e = #d c * c (4.23b) 
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in terms of the gravito-electromagnetic potential. The mass-current can be similarly split into its timelike 
and spacelike parts with respect to £ as 



where i^J^ = . 

The transverse gauge condition in terms of £, ip£, and is: 

V-V = S^ b + do. - V^j V-E, - B,r, = 0. 

For a gravito-electromagnetic frame £ (V£ = to order e) this reduces to 

V-V = S^C + d<^ - V ? ^ e - V-S c = 0. 
Using d4.22|) , flA.19p and ( |A.3C| ) it follows that d^ - V^S = - = hence 

V-V>= (<5^^ b ) + (^-V-S ? ) = 0. (4.27) 
The two bracketed terms are orthogonal so the transverse gauge condition is equivalent to the two equations 



(4.24) 

(4.25) 
(4.26) 



Hi = o, 

£ € = V-S 5 . 



(4.28a) 



(4.28b) 

The first condition is analogous to the Lorenz gauge in electromagnetism, while the second condition has no 
electromagnetic analogue. The consequences of this second condition are explored below. 

The equation (4.28a) implies (Appendix |X]) that ( 4.1 2| ) takes the Maxwell-like covariant form 

6T S = -2kJ(. (4.29) 

The perturbative part of the stress energy-momentum tensor may be expressed in terms of the mass- 
current using (4.13), (4.15) and (4. IS) as 

f = P^ b ®^- Jz®?-?® J( + t 6 , (4.30) 

and since the background source is assumed zero (T = 0) the divergence condition becomes 

V-t = (4.31) 

hence 

(<5J?£ b )-(V^-V-T e )=0. (4.32) 
Since the two bracketed terms are orthogonal one has: 

6Jt = 0, (4.33a) 

V 5 J c = V-T e . (4.33b) 

with the first condition expressing the conservation of mass-current in the background geometry. 

The field equation ( 4.29 ) and the closure of ( 4.20 ) can now be written in terms of these decompositions: 



d e #B e = 0, 
d e £ ? + £ e #B 5 = 0, 
d e #^ = -2«^#1, 
d e B e - CtfSi = -2k#J ? . 



(4.34a) 
(4.34b) 
(4.34c) 
(4.34d) 



Although the equations (4.34) involving the gravitational field have the mathematical structure of Maxwell's 
equations for electromagnetism the two theories are not isomorphic since the gravitational fields must ad- 
ditionally satisfy (4.28b) in the transverse gauge. Not all solutions of Maxwell's equations translate to 
perturbative gravitational fields that are compatible with this condition for a given S^. There does exist 
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however a class of compatible solutions. Such solutions ip are characterized by the existence of a gauge in 
which 

V> = 0, (4.35a) 

= 0, (4.35b) 
in some gravito-electromagnetic frame £ (Section |4.1|). In this case such solutions will be said to belong to 



the gravito-electromagnetic limit. It follows from equation (4.1) that 

S« = &gf. (4.36) 

Such a £ will not be unique. If C = £ + v, such that <?(£, £) = — 1 and w 2 = <?(v, v) <C 1, then 11^0 is of order 
ve. Hence if v < e, C a l so defines the gravito-electromagnetic limit. However the fields (V>£, S^, &q ) 
defined with respect to £ are identical (within this approximation) to the corresponding fields defined with 
respect to £ (App endix |Q ) . 

^3b[) a 



The condition fl4.35b| ) and the field equation ( |4.10D require 

r £ = o. 

Equation ( 4.35b ) implies V-Sj = d^tp^, hence in this limit condition ( 4.28b| ) becomes 

£ 4 = d^, 



or with (|4.23aD 



(4.37) 
(4.38) 
(4.39) 
(4.40) 



and using ( 4.28a ) it follows that 

= o. 

Thus is time independent and it follows from ( 4.23b| ) that the gravito- magnetic field is also time inde- 
pendent, C(15^ — . In this limit the field equations reduce to 

d 4 #B c = 0, 



d^ = 0, 
d € B € - £ s #£ s = -2k#J 5 , 



along with the conditions ( 4.28a ), ( 4.28b ) which can be written as 

d ? #*? - = 0, 

£ ? * ? = 0. 



respectively. Similarly it follows from ( 4.33a ), ( 4.33b ) and ( 4.37 ) that 

d ? #J ? +£p C -0, 



0, 



(4.41a) 
(4.41b) 
(4.41c) 
(4.41d) 

(4.42a) 
(4.42b) 

(4.43a) 
(4.43b) 



With £ = and the vector fields E^ = £|, = B^, JJ^ = J^, the field equations ( 4.41 ) can be written 
in more familiar notation as 



divl 



0, 



curlEfr = 0, 
divEj = —2kp^, 



curl I 



<9E 5 
~dt 



-2k% 



(4.44a) 
(4.44b) 
(4.44c) 

(4.44d) 
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dt 



divAj — 
9Aj- _ 

and the conservation equations ( 4.43a| ), ( 4.43b| ) become 



divJe + ^ = 0, 



dt 



0. 



(4.45a) 
(4.45b) 

(4.46a) 
(4.46b) 



Thus stationary electric and magnetic field configurations and their sources have direct analogues in the 
theory of perturbative gravitation. 

The Newtonian limit is defined as the gravito-electromagnetic limit supplemented by the condition 



From ( 4.12 ) this implies 

The Newtonian potential <£>£ is identified as 



4V 



0. 



0. 



From (4.42a) it must be time-independent 



and equation fl4.14|) becomes 



A$£ = 47rG/9£, 



(4.47) 
(4.48) 
(4.49) 
(4.50) 
(4.51) 



which is just the field equation for Newtonian gravitation. 

In the following section the above framework is illustrated by calculating the gravito-electromagnetic fields 
arising as perturbations on the asymptotic gravitational field of a rotating source. 

4.3. The gravito-electromagnetic Field of a Rotating Source. 

The metric tensor at large distances from a compact rotating body, with Newtonian gravitational mass 
M and angular momentum J, may be approximated (for r 3> 2GM in units where c = 1) by 



<J 



1 - 
2GJ 



2GM \ 



dt <g> dt 



1 



2GM \ 



• sin fir sin 



r J 

< dt + r sin 9dt <g> d(j>). 



(dr (g)dr + r 2 d9(g>d9 + r 2 sin 2 



(4.52) 



In the cobasis {e° = dt, e 1 = dr, e 2 = rd9, e 3 = r sin 9d<j)} the metric tensor can be rewritten as 



= -u-*h. 

r 



o 



1 + t)< s ' 



e 1 + e 2 ® e 2 + e 3 < 



e 3 ) 



^isin0( e o ®e 3 + e 3 ®e ) 



(4.53) 



in terms of the two length-scales i?s = 2GM and i?K = jj ■ 
Defining the dimensionlcss coordinates R and T by 

e 



Rs 



T, 



(4.54a) 
(4.54b) 



14 



S. J. CLARK AND R. W. TUCKER 



and the dimensionless cobasis e a by 



the dimensionless physical metric tensor g defined by 



Rs 



2 



9, (4.56) 



can be rewritten as 

9= - (l - ~J g° <8> e° + (l + ^) (e 1 <g> e 1 + e 2 <g> e 2 + e 3 <g> e 3 ) 



+ -^Asin6>(e° <g e 3 + e 3 <g> e°), 



(4.57) 



where the constant A is 



A = e|i. (4.58) 

For the last term in ( 4.57 ) not to be second order in e (and hence negligible), j^ 3 - must be of order e. 
With the dimcnsionlcss Minkowski metric tensor 

g = -e° <g> e° + e 1 <g> e 1 + e 2 ® e 2 + e 3 <g> e 3 (4.59) 

the perturbation /i about g is then 

ft = — (e° «) e° + e 1 (g) e 1 + e 2 ® e 2 + e 3 ® e 3 ) + A^A sin6>(e 3 <E> e° + e° ® e 3 ), (4.60) 
it it 2 

and taking the trace-reverse yields 

ip = 2-ie <g> e° + A sin (9(e 3 ® e° + e° <g> e 3 ). (4.61) 
— R R 

Using the dimensionless basis {X Q } dual to {e a } given by {2£o = dr,2Li — dn,2L2 — "^5e,X 3 = Rf \ e d^}, 
we use 2L.Q to define the dimensionless gravito-electromagnetic 1-form potential tp as 

i> x =2^e° + -^Asin0e 3 . (4.62) 
— A K R z 

In terms of the dimensionless gravito-Faraday 2-form £jc defined by 

the dimensionless gravito-electric and gravito-magnetic fields, £_ x an d follow by writing 

E 2Lo =^hS Kjy +±B Kn , (4.64) 



Hence from ( 4.62 ) 

and (with #1 = e 1 A e 2 A e 3 ) the dimensionless gravito-magnetic field is 



^ = -2^, (4.65) 



% = ^3 A ( sin61 ^ 2 + 2 c os Be 1 ). (4.66) 
In terms of Z = i?cos#, | the metric-dual map associated with g and 

A = Ad z , (4.67) 

the gravito-electric and -magnetic vector fields defined by M.x = £-x an d Sjc = become 

%„ = - 2 ^i' ( 4 ' 68a ) 

%„ = -^3 ( A " %( A ^i)^i)- ( 4 - 68b ) 

These fields may be compared with those derived in pj| and |24]. In the next two sections the motion of 

a test particle in a gravito-electromagnetic field is discussed. For such a particle with "spin" the motion 
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is expected to follow from the weak field limit of equations presented in p). To facilitate a comparison 
with treatments that ignore Mathisson-Papapetrou type coupling to spacetime curvature the discussion is 
restricted to the separate pre-geodesic motion of a spinless test particle and the gyroscopic precession that 
follows from a parallel spin vector along such a motion. 

5. Motion of a Test Particle in a Weak Gravitational Field 



The history of an electrically neutral, spinless test particle in spacetime is modeled by a future-pointing 
timelike curve C(t) that satisfies 



(5.1) 



2 g{C',C) 

for some general parameter r. To compare the following gravito-electromagnetic equations with analogous 
equations from electromagnetism in Minkowski spacetime it is most convenient to parameterize the curve so 
that 



g(C',C) 



-1. 



Using (3/7) and ( 3.11 ), to first order in e equation ( |5.l|) is 



V C >C + j(C', C, -) - -C'{h{C, C'))C = 0. 



(5.2) 

(5.3) 

(5.4) 
(5.5) 



Rewriting 7(6", C", — ) in terms of Tc and £c" (see Appendix |B|), (|5.3| ) becomes 

Vo'C" + (ic^c + \&C> Tr(S c 0) S = 
and, for a test particle with mass m, the perturbed gravitational force Fq, is then 

F G = -m(i c >T c > + |dc Tr(S c 0)'. 

The first term on the right is analogous to the Lorentz force in electromagnetism, however the second term is 
non-Maxwellian. Since Fq depends only on Tc and He it is invariant under gravito-electromagnetic gauge 
transformations (Section p~l] ). 

To examine the above equations in the gravito-electromagnetic limit assume that some vector dt defines 
a gravito-electromagnetic frame, and reparameterize C(r) in terms of t, so that C(t) = C(r). Writing C'(t) 
as 

C' = d t +v, (5.6) 

where g(v, dt) — 0, it then follows that 

g(C',C')=-(l-v 2 ), (5.7) 

where v 2 = g(v,v). 
Writing C 

so that v 1 (hence terms smaller than ev will be neglected) , ( |5.4|) can be written in terms of C as 



-^==yC', and introducing the assumption that the speed v of the particle is non-relativistic, 



VcC + vd t vC + (icFa + jd a Tr(S c ,)) J = 0. 
Expressing Tc and Sc/ in terms of ipg t (see Appendix ^ with C = C and ^ = dt) so that 

T C '=T dt +d dt (i> dt (v))Adt, 



it follows that 



Sc = {<fe t + 2^a t (v)}ga t - ip dt ®v b + v b ® ^ 9t , 



ie>T c > = -8* - #(v b A B a J - d 9l (^ (v)) + £ 9l (v)dt, 
dc- Tr(S c = 3£ dt + 4d at (^8»(v)) + 33 t 9t v - 3£ 9t (v)di, 



(5.8) 

(5.9a) 
(5.9b) 

(5.10a) 
(5.10b) 



where (1.38) has been used. 
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— = 3^+#KAB 8t ) - 7 ^ at v, (5.11) 



Using (5.6) the spatial part of the equation of motion (5.8) becomes 

^-K+^a*)'-! 

and if $>a t is time-independent this reduces to 

dv 



or in vector notation 



^-l4 ( +#( v AB 0t J . (5.12) 
dv 

— = iEa t +vxB St , (5.13) 

where Eg, = £ a and Ba f = &g ■ 

In terms of these fields the perturbed gravitational force Wq is then 

F G = ™(±Ea t +vxBa t ). (5.14) 

This is similar in structure to the Lorentz force law of electromagnetism, except for the factor of j multiplying 
the gravito-electric term (and the fact that the Ea t and Ba t fields couple universally to inertial mass). Note 
that rescaling Ea t to remove the factor j would introduce a factor int o the g ravit o-Maxwell equations above. 



If one continues to the Newtonian limit then Ba t = and, using ( 4.38 ), ( 4.4S ) to express Ea t in terms of 



the Newtonian potential $a t (which will automatically be time- independent by the gauge condition V-?/? 



0), fl5.14|) reduces to 

F G = -mgrad$a t , (5.15) 



which is just Newton's force of gravity in terms of a gravitational potential <!>at satisfying (4.51) 



6. Precession of a Small Gyroscope in a Weak Gravitational Field 
One may model the relativistic spin of a freely falling gyroscope by a unit spacelike vector field S along 



a future-pointing timelike curve C(r) that satisfies (5.1), such that S solves 



V c <S = 0, (6.1) 

with 

$(S,S) = 1, (6.2a) 

ff(S,C") = 0, (6.2b) 

in terms of the physical metric tensor. A perturbative analysis can be given in terms of the vector s defined 
along C(t) by 

s = (1 + I/i(S,S))S-/i(S,C")C", (6.3) 

so that (to first order in e) 

ff(s,s) = l, (6.4a) 
g(B,C)=0, (6.4b) 

in terms of the background metric g. 

Equation ( |6.3[ ) can be inverted and to first order in e 

S = (1- ±h(s,s))s + h(s,C')C (6.5) 

Using the perturbed connection, fl6.l[) can be written in terms of s to first order in e as 



V C 'S - ~C'(/i(s, s))s + C'(h(s, C'))C + 7(C", s, -) = 0. (6.6) 
Rewriting 7(6", s, — ) in terms of Ta and Sc (see Appendix |b|), this becomes 

v£,s + i(n c d s jF c , -n s v c -{s c -( s ,-)})» = 0, (6.7) 



GAUGE SYMMETRY AND GRAVITO-ELECTROMAGNETISM 



17 



where for any vector field X, the Fermi- Walker connection V F is defined on C by 

V F C ,X = V C >X + g{C, X)Ac> - g(Ac- , X)C, 



(6.8) 



(6.9) 



with A-c = VcC". From (6/7) the effective gravitational torque Tq on the gyroscope is 

T G = -~(IL c >i s Fc> - n s V C '{£c<(s, -)})'. 

The vector s is said to be non-rotating along the path C when Tq = 0. As with the point particle, if 
this curv e is parameterized in terms of t so that C'(t) is given by (5.6) (and dropping terms smaller than 
ev), (6.7) becomes 



v£, s + i(n c ,i s jr c , - n s v c ,{£ c ,(s, -)})« = o. 

The vector s can be written in terms of its spacelike component <r as 

s = g((T,v)d t + <x, 



(6.10) 



(6.11) 



which follows from (3.4a). Imposing the gravito-electromagnetic limit in the dt frame, and writing Tc and 
Sc in terms of Eg t , Mg t and <f>g t (see Appendix |c|) , the spacelike component of ( 6.1C ) then becomes 



der 

dF 



i (tr x B ft + (a- • v)E 9t - i(o- • E a Jv + (C> 9 > 



(6.12) 



Since the length ^jcr ■ er of the vector a is dependent on t (to the appropriate order), introduce the 
constant length gyroscopic spin vector S by 

$ = (l-i<£a>-!(v-o-)v. (6.13) 



In terms of S, (6.12) can be written 



— = -S x 
di 2 



- -v x Ea t 



(6.14) 



To this approximation the precession rate 5 (§v x Eg, — Bg t j of § is independent of the gyroscopic spin 



which may be compared with the similar equation found in jllj 

To this approximation the precession 
and could in principle be used to detect Mg t . 

7. Field Redefinitions 



In order to facilitate a comparison with alternative formulations of weak gravity using the gravito-electro- 
magnetic analogy it is worthwhile to effect certain field redfinitions. These are motivated by looking at the 
analogues of the transformations (3.37) with (4.5). 

Let = /i(£, — ) and £ be a gravito-electromagnetic frame, then under a gravito-electromagnetic gauge 
transformation transforms as 



h(. h-> ^ + dA-(£A)f , 



which is not quite of the same form as (2.2). However, defining 



and 



it follows, writing £ = dt, that 



(7.1) 

(7.2a) 

(7.2b) 

(7.3a) 
(7.3b) 



A h + grad A 

which is indeed analogous to the (3+l)-decomposed form of an electromagnetic gauge transformation. 
Similarly ip£ transforms as in ( |4.7[ ) and defining 

** = -& = -lK£,0, (7.4a) 
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the (3+l)-decomposed form of the gauge transformation ( |4.7| ) becomes 



L </> 



grad A. 



(7.4b) 

(7.5a) 
(7.5b) 



Since (7.3), (7.5) are analagous to electromagnetic gauge transformations this suggests that a Maxwell 
type syste m o f field equations may be constructed using <E>/j and A/j instead of <£",/, and A^,. 
^From (11) and the above definitions of Ah and A^, it follows that 



$V, = $ ft -lTr(S), 
A v , = A h , 



(7.6a) 
(7.6b) 



where £ = S^. 

The Lorcnz gauge condition in ( 4.28a ) written in terms of $rt and Art becomes 

divArt + d t $v; = 0, (7.7a) 

which has the same form as in electromagnetism for all S. In terms of and Ah this becomes 

dWA h + d t <P h = ±Tr(£), (7.7b) 

which is unlike the electromagnetic Lorenz gauge condition when the right-hand side is non-zero. 

The gravito-electric and -magnetic fields E,/, = and = discussed earlier are related to <!>,/, and 
Ay, by 

Ert = -grad $rt - d t A^, 



Brt = curl A^. 

Alternative gravito-electromagnetic fields Eh and can be defined in terms of $>h and Ah as 

E/j = -grad$ ft - d t A h , 



These are related to E^ and B^ by 



E,j 



B, 



B, 



curl A/, 



igradTr(S), 



B, 



In terms of E^ and £ the second gauge condition ( 4.28b| ) can be written as 



(7.8a) 
(7.8b) 

(7.9a) 
(7.9b) 

(7.10a) 
(7.10b) 

(7.11a) 

(7.11b) 

(7.12a) 
(7.12b) 
(7.12c) 

>rt - dtErt = -16ttGJ (7.12d) 

for all S. As stressed above these equations together with ( 1.17| ), ( 7.7a ), ( [7. 11a ) summarise linearised 
Einsteinnian gravitation. 



or using E^, and £ as 



Ert = (V-£) B , 
E, l = (V-S) S -igradTr(S). 



In the E^ and B^ notation , with k = 8irG, p = p^ and J = 

div B^ = 0, 

curl E^ + <9*Brt = 0, 



the field equations (4.34) become 



div E 



rt 



-\QirGp, 



curl B„ 



Alternatively using ( 7.10a ) and ( 7.10b ) the above take the form 

divB h = 0, 



(7.13a) 
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CUrlEfc + dtBft = 0, 

divE/, = -WnGp- ±divgrad Tr(£), 

curlB/, - ftE h = -16ttGJ + ±<9 t grad Tr(£) 

in terms of E/, and B/j and S. In the gravit o-ele ctromagnetic limit, Tr(£) 
E„ 



-12% so $,/, 



4E/j. Then the Lorenz gauge condition (7.7b) becomes 



divAft + 4d t <S>h = 0, 



and (since c^B^ = in the gravito-electromagnetic limit) the field equations ( 7.13 ) simplify to 

divBft = 0, 

curlE/j = 0, 

div E h = -AixGp, 

curlB h - 4<9 t E /l = -16ttGJ 
which may be compared with (4.44) and the formulations given in f23|| , pl| . 

8. Discussion 



(7.13b) 

(7.13c) 

(7.13d) 
4$/j and 

(7.14) 

(7.15a) 
(7.15b) 
(7.15c) 
(7.15d) 



In this article the analogy between Maxwell's equations for the electromagnetic field, (2.1a), (2.1t) and 
the Einstein equations for weak gravitational fields in the transverse gauge, (4.20), (4.29) has been made 
in terms of tensor fields. While the former are valid in an arbitrary Lorentzian spacetime the latter have 
been developed in terms of perturbations about a flat spacetime background. A comparison has been made 
between the general equation describing the motion of electrically charged point particles (2.8) and the motion 
of massive point particles in a weak gravitational field (5.4). Equations have also been developed (6.14) 
describing the motion of a freely falling small gyroscope in terms of gravito-electromagnetic fields. 

In general it is asserted that any analogy between electromagnetism and weak gravity is closest in a 
restricted class of reference frames related by suitable non-relativistic transformations and for stationary 
physical field configurations in such frames. In addition to gravito-electromagnetic fields the general equations 
of weak gravity involve a second degree symmetric tensor field which has no electromagnetic analogue. 
The gravito-electromagnetic fields defined in 4.2 are coupled to via ( 4.28b| ) and this field produces non- 
Maxwellian terms in the weak gravitational force and torque equations, (5.5) and fl6.9|) respectively. 



In electromagnetism the Maxwell fields A, F and J are defined independent of any frame of reference. The 
latter is only re quired t o def ine electric and magnetic fields and their sources in terms of electric charge and 
current density ( 2.10 ), ( 2.11 ). In gravito-electromagnetism the analogous fields ipg, and are manifestly 
f r ame- dependent. 

The definition of these fields has been motivated by their behaviour under a class of gauge transforma- 
tions belonging to the gauge symmetry of the weak Einstein equations. Unlike electromagnetism, is not 
gauge invariant under these transformations in general. However a subset of these transformations does exist 
for which J-^ remains invariant. These have been called gravito-electromagnetic gauge transformations by 
analogy with the gauge symmetry of Maxwell's equations. Unlike electromagnetic interactions with elec- 
trically charged particles, weak gravitational interactions are not mediated by complex representations of 
these symmetries. In terms of the gravito-electromagnetic fields a subset of the linearised Einstein equa- 
tions take a remarkable form that is isomorphic to Maxwell's electromagnetic field equations. No limits or 
further approximations are required to establish this correspondence. The explicit appearance of the 
tensor is confined to the remaining equations in the linearised system. This is a primary distinction of the 
approach adopted here compared with previous derivations of the gravito-electromagnetic field equations. 
To exploit this reformulation and link physical field configurations with solutions to Maxwell's equations 
further conditions must be imposed on the linearised system. 

Conditions have been found that enable a useful analogy between weak gravitation and electromagnetism 
to be established. In the gravito-electromagnetic limit depends only on the a gravito-electromagnetic 
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potential ( |4.36| ). Consequently, such weak gravitational fields can be described in terms of ip^ (or alternatively 
and B$). 



However as stressed above, in order to obtain Maxwell-like equations the transverse gauge condition (4.11) 
is imposed on ip. Condition (4.27) induces an equivalent condition ( 4.28a| ) on ip^. (By contrast the electro- 
magnetic gauge condition (2.6) is one of many that may be imposed on A.) The condition on ip also imposes 
the restriction (4.28b) which in the gravito-electromagnctic limit implies that ip£ has a restricted time depen- 
dence. Physical gravito-electromagnetism consequently shares more in common with electromagneto-statics 
than electromagnetism. Although both the Maxwell equations (2.12) and the equivalent gravito-electro- 



magnetic equations (4.34) hold in any inertial frame (in Minkowski spacetime), in order to remain within 
the gravito-electromagnetic limit only a class of gravito-electromagnetic frames of reference is permitted 
(Appendix ||) . 

By perturbing the equation of a physical timelike geodesic the relativistic equation of motion for a massive 
point particle in t he weak gravitational field can be cast into a form containing a gravito-electromagnetic 
Lorentz force ( |5.5| ) and an additional non-Maxwellian term. It is worth pointing out that in the context of 
the gravito-electromagnetic fields defined in this paper the derivation of this equation of motion does not 
rely on the speed of the particle and the gravito-electromagnetic Lorentz-like force takes its natural form. 
In the gravito-electromagnetic and non-relativistic limit the particle acceleration is then determined by a 
non-relativisti Lorentz-like force ( 5.14 ) containing an additional factor of \ multiplying the gravito-electric 
field. 

When working in the gravito-electromagnetic limit the field redefinitions presented in section |t] permit a 
comparison with the work of Thorne in [^3| with the notation $/j = $, Ah = 7, E;,, = g, and B^, = H, and 
with a current of the form J = pv. The work of JlTQ , fll6| , and others may be related to that of |2j| 

either by trivial field redefinitions or changes in metric signature. The analogy between electromagnetism and 
weak gravity developed by Wald |2[| is similar to the one presented here. However he does not discuss how 
the gravito-electromagnetic 1-form-potential behaves under gauge transformations nor how the restricted 
time dependence arises from the transverse gauge condition. 

This article offers an alternative description of weak field gravitation in the language of gravito-electro- 
magnetic fields. We feel that it clarifies a number of issues concerning various other analogies between the 
equations of post-Newtonian gravity according to Einstein and Maxwell's description of electromagnetism. 
In the absence of a gravito-electromagnetic limit formulations based on <&h and give rise to gauge 
conditions ( 7.7b[ ), ( 7.11b| ) and field equations (7.13c), (7.13d) thereby exposing couplings between the gravito- 
electromagnetic potentials and fi elds an d the tensor field . In the approach adopted here such couplings are 
relegated to the gauge condition ( 7.11a ). Many analogies coalesce in the gravito-electromagnetic limit modulo 
a re-shuffling of numerical factors that cannot be scaled away entirely. Although the mathematical analogy 
between weak gravity and the full system of Maxwells equations for electromagnetism in terms of covariant 
tensor fields on flat spacetime can be made close, the existence of gravitational gauge conditions limits the 
physical^ analogy to stationary phenomena. Despite this limitation the interpretation of weak gravity in 
terms of gravito-electromagnetic fields offers a fertile avenue of exploration for phenomena associated with 
the detection of the stationary gravito-magnetic field. The methods presented here are also applicable in 
principle to certain non-flat backgrounds and to weak field descriptions of non-Einsteinian gravitation (in 
which other geometrical fields may compete with metric-induced gravity) at the post-Newtonian level. These 
issues will be discussed elsewhere. 
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Appendix A. Definitions and Notation 

Throughout this article the geometry of spacetime is described in terms of a smooth Lorentzian metric 
tensor field and its associated Levi-Civita connection. The connection 1-forms ui a b associated with any such 
connection V are defined with respect to any local basis of vector fields {X a } where a = 0, 1, 2, 3 by 

V Xa X b = u c b {X a )X c . (A.l) 

If T is a smooth tensor field on spacetime then in terms of this connection the covariant differential V of 
T is defined by 

(VT)(A, -,...,-) = (V*T), (A.2) 

and if the first argument of T is contravariant, the divergence of T is 

\7-T=(X7 x T)(X a , -,...,-), (A3) 

where indices are raised with the components of the metric tensor in the usual way. Similarly with {e a } dual 
to {X a } (i.e. e a (X b ) = 5%) 

VS=(V Xa S)(e a , -,...,-), (A.4) 

for all tensors S whose first argument is covariant. For any degree 2 covariant tensor T the tensor SymT is 
defined by SymT{X,Y) = \{T{X, Y) + T(Y, X)). 

In terms of the exterior derivative d on smooth differential p-forms on spacetime, the covariant exterior 
derivative is defined on mixed basis indexed p-forms by 

D5°- 6 c ...d =dS"*- 6 c ...d + Lu a s A S'- b c.. d + ... + Lu b s A S a - S c ... d 

a b sab ( A ' 5 ) 

— lo c A S "' s ...d — ... u) jAS "' c ... s . 
The curvature operator R for the connection V is 

Rxy = [V x ,V y ] - V [X ,Y] (A.6) 
for all vector fields X, Y, and the (Riemann) curvature tensor R defined by 

R(X,Y,Z,a) = a(R X YZ). (A.7) 
Then in dual local bases {X a } and {e a } the curvature 2-forms R a b are defined by 

R = 2R a b ®e b ®X a . (A.8) 
Successive contractions (where i x is the interior operator with respect to the vector X) give the Ricci 1-forms 

P b = ix a R\, (A.9) 

and the curvature scalar 1Z 

Tl = l X aP a . (A.10) 

With g a metric tensor on vectors and G the induced metric tensor on 1-forms, the metric-dual of a vector 
X is given by X b = g(X, — ) and that of and a 1-form a by a" = G(a, — ). 

If U is a unit-normalized vector field (g(U,U) = A where A = ±1), then the projection operator on 
contravariant tensor fields is defined as 

Uu = 1 - \U®U b . (A.lla) 
By abuse of notation the same symbol is used to denote the projector on covariant tensors 

Tlu = 1 - AC/ b <g> U, (A.llb) 

and differential p-forms 

lit/ = 1 - XU b Ahy, (A. 11c) 
since the domain should be clear from the context. The map Tlu is a tensor homomorphism 

n I7 (a<g>...<g>/3) = n I 7a<g>...<g>n { 7/3, (A.lld) 

for all a, (3, ... and has the following properties: 

UuUu = Uu, (A.12a) 
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U u i u = i u , (A.12b) 

iuUu = 0. (A.12c) 

The projection tensor Ily is used to define the spatial metric on spacetime associated with the timclikc 
vector field V: 

gv = U v g. (A.13) 
One can write W b in terms of gy and Ily as follows: 

W b = (t v + fl v + ±6ygy - v b <g> Ay ( A - 14 ) 

where 

try = SymIVW b -±e v gy, (A.15a) 

is the shear of V, 

9y = V-V, (A. 15b) 

is the expansion of V, 

fl v =ttvdV b , (A.15c) 

is the vorticity of V" and 

Ay = VyF (A.15d) 

is the acceleration of V. 

A vector field £ is said to be parallel (with respect to V) if 

V£ = 0, (A.16) 

in which case V£ b = 0, so £ has vanishing shear, vorticity, expansion and acceleration. Furthermore 

SymV£ b = 0, (A.17) 

or equivalently 

= (A.18) 
in terms of the Lie derivative. Thus, if £ is parallel then it is also a Killing vector and 

V € = £ € . (A.19) 
For any tensor field T taking at least one covariant argument let be defined by contraction such that 

7H-, ...,-)= T(- ...,-,£,-,...,-). (A.20) 

If £ is parallel it follows that 

(V*T € )(- ...,-) = (V*T)(- ...,-£,-...,-) (A.21) 

for all vector fields X . When acting on p-forms this is simply the rule: 

i«V x = V x i 5 . (A.22) 

The metric tensor g gives rise to a canonical volume 4- form *1 on spacetime and an associated Hodge map 
* on p-forms. In terms of a local g— orthonormal local basis of 1-forms {e a } one may write *1 = e°Ae 1 Ae 2 Ae 3 . 
A volume 3-form #1 associated with the unit timelike vector field V is given in terms of *1 by 

*1 = V b A #1 (A.23) 

with iy#l = 0. Hence 

#1 = - * V b (A.24) 

and #1 induces a spatial Hodge map, #, on the image of forms under Ily. With these operations defined 
one can decompose the Hodge map of any spacetime form and its exterior derivative into spatial forms. Thus 
the Hodge dual of any p-form u on spacetime may always be written 

*UJ= {#iyw} + V b A {#(nyw)"} (A.25) 
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where a n — (— l) p a for any p-form a and each of the terms in brackets is annihilated by V (i.e. iy| j = 0). 
A projected Lie derivative with respect to V is defined as 



and for any p-form oj one can write 

LyLO = |£yw — Ay A iyw j — V b A |£yiyw j 

where each of the terms in brackets is annihilated by V . 
In terms of the projected exterior derivative 



and with Dy defined by 



one may write 



dy = Ilydlly, 



Dy = dy + Ay A 



duj = |dyCL> — f2y A iyw I — V b A |(£yCJ — Dyiyw) j 



(A.26) 



(A.27) 



(A.28) 



(A.29) 



(A.30) 



where each of the terms in brackets is annihilated by V. These formulae permit a local "3+1" decomposition 
of exterior differential equations with respect to the general observer field V in a spacetime with an arbitrary 
Lorentzian metric and permit one to identify spatial fields parametrised by a local time associated with V . 
With the Faraday 2-form F decomposed as 



F = V b A e + #b 

where iye = 0, and iyb = 0, it follows that iyF = — e and TlyF = #b. Similarly using ( A.25| ) 

*F = V b Ab- #e, 

so iy * F — — b and IlyF = #e. Writing 

J = pV b +j 



where iyj = 0, and using (A. 25) it follows that 



*J = -P#l - V b A #j. 



(A.31) 



(A.32) 



(A.33) 



(A.34) 



The co-derivative 8 is defined on spacetime p- forms in terms of the Hodge map * and the exterior derivative 
dby 



S = * 1 d * if 



where 7]uj = to 71 for any p-form lj. The field equations (2.1a) and (2.1t) can be written as 

dF = 0, 



(A.35) 



(A.36a) 



d* F = *J. 

Using (A.3C) the Maxwell equations given in (2.12) follow immediately. 



(A.36b) 
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With the definitions: 

divX= (#d v #X b ) J , (A.37a) 

curlX= (#d v X b ) S , (A.37b) 

grad$ = (d y $) S , (A.37c) 

XxY=(x b AY b ), (A.37d) 

X-Y = g y (X,Y) = #(X b A#Y b ) (A.37e) 

where X and Y are spacelike (with respect to V) vectors, and $ is a 0-form on spacetime, exterior equations 
can be transcribed to Euclidean vector notation. 

The Laplacian operator Lap associated with g and V is defined by 

Lap = V-V, (A.38) 

and the associated Laplace-Beltrami operator A on p-forms is 

A = -(Sd + dS). (A.39) 

It can be shown that, when acting on any p-form a, Lap is related to the Laplace-Beltrami operator A by 

Aa = Lapa + e a Ai X bIlx a x b a. (A. 40) 

Thus if (3 is a 1-form 
and for any 0-form / 
In a spacetime with a flat metric 
for a/Zp-forms. 



A/3 = Lap/3- Ric(/3 B ,-) (A.41) 
A/ = Lap/. (A.42) 
A = Lap (A.43) 
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Appendix B. The Perturbed Connection and Gravito-Electromagnetic Fields 
With the (3,0) tensor field 7 b defined by 

7 b (x,y,z) = 7 (x,y,z b ), (b.i) 



using ( 3.17 ) it follows that 

7 b (C, -, -) = §V e /i + Alt(Vft)(- f , -), (B.2) 

where £ defines a gravito-electromagnetic frame (unit timelikc and parallel, to zero order in e), and for any 
type (2,0) tensor T, Alt T is defined by Alt T(X,Y) = \{T{X,Y) - T(Y,X)). 

To relate ex pres sion (B.2) to gravito-electromagnetic fields one may proceed as follows. By trace-reversing 
ip as given in (4.1), h can be written as 

h = ® £ b - C b ® ^ ~ S e + i Tr(S ? ) ff , (B.3) 

and acting with gives 

V $ h= - ® £ b - £ b <8> Vf^£ - V 5 S 4 + ±£ Tr(S e )sf. (B.4) 
Since V£ is first order in e it follows that 

(V/iX- £,-) = V(/i(£,-)), (B.5) 



(B.9) 



and contracting (B.3) with £ gives 

=^-(^-!Tr(£ ? ))£ b . (B.6) 

Hence 

(V/i)(-, £, ") = - d(fc - \ Tr(S c )) ® £ b . (B.7) 

Antisymmetrizing this yields 

Alt(V/i)(- £ -) = ^ + e b A d 4 (^ - i Tr(S c )), (B.8) 

where = d^. 

Equation (B.2) can now be written as 

7 b (C, -, -) = - (v^ e ® t + £ b ® v^j + v € s s ) 

+ ^Tr(S ? ) 5 + ^ + C b A - § Tr(S ? )), 

and contracting on £ gives 

7 b (£, t -) = -e(0« - | Tr(S 5 ))^ + + id c Tr(S ? ). (B.10) 
If X is orthogonal to £, <?(£, X) = 0, and VX is first order or higher in e, then 7 (£, X, — ) can be rewritten 

as 

7 b (C, X, -) = - £(^(X))£ b + £(i Tr(S c ) + E £ (X, X))X b ^ 

+ i x (i^ + |d ? Tr(£ e ))£ + |n s i x J^ + n x V e S ? (X, -). 

Let £ be a vector field such that V£ = and g (£,,£,) = —1 (to at least first order in e). Using (3.11) and 
adding hats to quantities defined with respect to the physical metric g it follows that 

V£ = 7(e."-). (B.12) 

or with X' b = G(X,-), that: 

V?W(£,-,-), (B.13) 
where 7 b (A", Y, Z) = 7 (X, Y, Z v ). Antisymmetrizing and using ( |B.S| ) yields 

d£ b = ^+£ b Ad c (^-±Tr(£ ? )). (B.14) 
The left hand side can be contracted on £ and rewritten as 

\^i = A £-£d(Mf,0), (B.15) 
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where A. £ = V^. Contracting ( |B.6| ) with £ and using ( B.14[ ) it follows that 

A ^-£ 5 + ±d c Tr(£ c )) + ±£Tr(£^. 



Defining 



and projecting ( B.14 ) gives 



« 4 = #B«- 



In the gravito-electromagnetic limit and with = ( B.lq ) simplifies to 



(B.16) 
(B.17) 
(B.18) 

(B.19) 



Thus the gravito-electromagnetic fields can be interpreted in terms of the vorticity and acceleration of 
the vector field £ with respect to g. 
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Appendix C. Transformation Laws For Gravito-Electromagnetic Fields Under a Change 

of Frame 



Let £ and £ define two gravito-electromagnetic frames (as defined in section 4.1) such that 

C = £ + v, (C.l) 

where v) = and v 2 = g(v, v) <C 1. 

Since tp is independent of any frame of reference it may be expanded as 

^ = <^S-^®£ b -£ b ®^-S c , (C2a) 

in terms of £ or as 

^ = ^0-^®^-^®^-^ (C.2b) 
in terms of £. Equating these two expressions yields 

A = ^ " ^c(v)£ b + {^v b - S c (v,~)} , (C.3a) 

S c = S ? + 2-0 5 (v)g c - ® vb - vb ® 

- {^v b - E c (v, -)} ® e b - ? b ® {0 ? v b - S e (v, -)} , 

Thus if v is of order e the fields ip^ and are invariant under a change of gravito-electromagnetic frame to 
first order in e. 

If the gravito-electromagnetic limit is satisfied in the £ frame, the terms in braces vanish and the above 
simplify to 

V>c=V^-V^(v)£ b , (C.4a) 

E c = + 2^(v)}g e - ipi ® v b - v b <g> V«- (C4b) 
If v is of order e then the gravito-electromagnetic limit is preserved and the above fields remain invariant. 



(C.3b) 
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Appendix D. An Alternative Analogy Between Einstein's Equations and 

Electromagnetism Based on Properties of the Conformal Tensor 

As mentioned in the introduction there exist alternative analogies between Einstein's theory of gravitation 
and Maxwell's electromagnetic equations that do not necessarily require a perturbative approach. One such 
analogy is summarised here in terms of the conformal tensor on spacetime since it highlights the differences 
accorded to gravito-electromagnetism by different authors [jl5| , , , || , || , Q . 

The fourth order (Weyl) conformal tensor may be written as 

C = —2C a b (8 (e a A e b ) (D.f) 

where (in any local cobasis {e a }), the conformal 2-forms C a b are defined in terms of the curvature 2-forms 
R a b, the Ricci 1-forms P a , |l| and the curvature scalar 1Z: 

C a b = Rab - ^{P a A e b - P 6 A e a ) + ^TZe a A e b . (D.2) 
The covariant exterior derivative of the conformal 2-forms is 

VC ab = ~(Y a Ae b -Y b Ae a ), (D.3) 

where the 2-forms Y a are defined by 

Y a = BP a - IdTZ A e a . (DA) 



If the geometry of spacetime is determined by Einstein's equations then these forms may be related to the 
stress energy-momentum tensor T by introducing the 3-forms r a such that 

T ={*- 1 r a )®e a . (D.5) 

Since the connection is torsion-free: 

Y a = B(P a - ift A e a ). (D.6) 
It follows immediately from Einstein's equations written in terms of the Ricci forms P a |ij that 

Ya = |D (*- x t - i(i X c *-! Tc )e a ) , (D.7) 

where k = 8nG. 

In terms of the covariant Lie derivative Lx : 

L x = Di x + ixD, (D.8) 

one may show that 

L X aC ab = ^Y b (D.9) 

and 

L x «Y a = 0. (D.10) 
These relations may be compared with the equations for the Faraday 2-form F: 

SF = J, (D.ll) 
from which there follows the conservation of electric current, 

SJ = 0. (D.12) 
In any local cobase F can be written in terms of its components as 

F = F ab e a Ae b , (D.13) 



(compare with (D.l)) and J as 

J = J a e a . (D.14) 
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In terms of the vector-valued 0-forms F a t and J a (D.ll) and (D.12) can be rewritten as 



Lx«F ab = -^J b , (D.15) 
L X 'J a = (D.16) 



which are analogous to ( D.9 ) and ( D.1C ) 



The use of local coframes is not mandatory to see this correspondence. A purely tensorial formulation 
follows by writing 

V-C=-F, (D.17) 

where 

Y = e a ®Y a . (D.18) 

This is the analogue of the Maxwell equation 

VF=~J. (D.19) 

Similarly one may write 

V-y = 0, (D.20) 
to compare it with the conservation of electric current written in the form 

V-J = 0. (D.21) 
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